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ABSTRACT 

Recent observations of several 7-ray burst (GRB) afterglows have shown evidence 
for a large amount of X-ray line emitting material, possibly arising from ionized iron. 
A significant detection of an X-ray spectral feature, such as that found in the Chan- 
dra observation of GRB 991216, may provide important constraints on the immediate 
environment of the burst and hence on progenitor models. The large Fe Ka equiv- 
alent widths inferred from the X-ray observations favor models in which the line is 
produced when the primary X-ray emission from the source strikes Thomson-thick 
material and Compton scatters into our line of sight. We present such reflection spec- 
tra here, computed in a fully self-consistent manner, and discuss the range of ionization 
parameters that may be relevant to different models of GRBs. We argue that the pres- 
ence of a strong hydrogen-like Ka line is unlikely, because Fe XXVI photons would be 
trapped resonantly and removed from the line core by Compton scattering. In contrast, 
a strong narrow emission line from He-like Fe XXV is prominent in the model spectra. 
We briefly discuss how these constraints may affect the line energy determination in 
GRB 991216. 

Subject headings: gamma rays: bursts — radiation mechanisms: non-thermal — line: 
formation 



1. Introduction 

The detection of spectral signatures associated with the environment of a 7-ray burst (GRB) 
would provide important clues about the triggering mechanism and the progenitor (Meszaros & 
Rees 1998; Lazzati, Campana & Ghisellini 1999; Bottcher 2000). Observations with Chandra, 
ASCA and BeppoSAX have provided tentative evidence for Fe Ka line and edge features in at least 
five bursts; GRB 970508 (Piro et al. 1999), GRB 970828 (Yoshida et al. 1999), GRB 991216 
(Piro et al. 2000) and GRB 000214 (Antonelli et al. 2000) all show an emission feature during 
the X-ray afterglow a few hours to a day after the burst event; while GRB 990705 (Amati et al. 
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2000) displays a prominent X-ray absorption feature during the burst itself. Although most of 
the line detections are only marginally significant and fail to distinguish between the various line 
excitation mechanisms \ GRB 991216 shows a 3.49±0.06 keV line at a moderate confidence level 
(~ 4a). This is consistent with emission from H-like Fe (Fe xxvi) at the redshift of the most 
distant absorption system along the line of sight at z=1.02 (Vreeswijk et al. 2000). A similar 
observation, but with higher signal-to-noise, may be able to distinguish between the various line 
emission mechanisms and lead to the correct determination of the line energy. 

The large equivalent widths (~ a few keV) inferred from the X-ray features favor models in 
which the line is produced by reflection (Vietri et al. 2000; Rees & Meszaros 2000), rather than 
transmission. Any detection of emission features in the afterglow spectra some hours after the burst 
(such as in GRB 991216) therefore imposes strong constraints on the location and geometry of the 
optically thick reflecting material. Observing an X-ray line at a time ?obs after the burst implies that 
the emitting material must be located within a distance ~ c?obs/(l +z) from the explosion site, thus 
strongly limiting the size of the remnant. However, this gas cannot be optically thick along the 
line of sight because this would smear out the short-time variability of the burst radiation. These 
conditions point towards a strongly anisotropic environment from which a GRB is seen only if we 
happen to observe the system through a line of sight with low optical depth (Bottcher 2000; Vietri 
et al. 2000; Lazzati et al. 2001). 

Two types of reflection models have been developed to explain the origin of the X-ray emis- 
sion features. The first invokes the interaction of the primary X-ray emission from the afterglow 
with an Fe-enriched, Thomson-thick, asymmetric remnant, which Compton scatters X-rays into 
our line of sight (Vietri et al. 2000; Bottcher & Fryer 2000). This scenario requires a mass 
> 0.06Mq of Fe at a distance about 1.5 light days, possibly due to a remnant of an explosive 
event or a supernova that occurred days or weeks prior to the GRB (see Vietri & Stella 1998; 
in contrast with MacFadyen & Woosley 1999 which favors a simultaneous supernova explosion). 
The other type of model involves a long-lived (> 1 day) magnetar or accreting black hole with 
a continuing but decaying outflow that interacts with the stellar envelope at distances less than a 
light-hour (Rees & Meszaros 2000; Meszaros & Rees 2001). In this case only a small mass of Fe 
is required, and can be readily produced by the star itself. 

Under both interpretations, it is likely that reflection takes place in highly ionized surfaces. 
This can lead to strong Comptonization of the emergent Fe line, and other absorption and emis- 
sion features. In this paper, we present and discuss detailed, self-consistent computations of the 
temperature and ionization structure of a uniform slab of gas ionized by the incident radiation of 



'The energy of the Fe emission hne shifts up in energy from 6.4 through 6.7 to 6.97 keV as Fe is ionized, but 
Doppler shifts may confuse precise estimates of the observed energy. 
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a GRB and of the resulting reflection spectra. Our analysis applies to an optically thick, homoge- 
neous medium, significantly extending previous analysis in the optically thin regime (e.g. Weth et 
al. 2000; Bottcher 2000), which is hard-pressed to explain the observed Fe-line feature in GRB 
991216 (Vietri et al. 2000). We estimate the strength of the Fe Ka line that each model produces. 
Finally, we discuss the implications of these results for current and future X-ray observations. 



2. X-ray illuminated slabs 

We employ the reflection code developed by Ross, Weaver & McCray (1978) and updated by 
Ross & Fabian (1993). We consider the illumination of the first 12 Thomson depths of a infinite, 
uniform slab of gas by radiation with a power-law spectrum of photon index T incident at an angle ^? 
to the normal (Ross, Fabian & Young 1999). The incident radiation is treated analytically in a 'one- 
stream' approximation, while the diffuse radiation that results from both the Compton scattering of 
the incident radiation and emission from the gas itself, is treated using the Fokker-Planck/diffusion 
method of Ross et al. (1978). Once thermal and ionization equilibrium in the slab is found, the 
reflection spectrum is computed. We assume that hydrogen and helium are completely ionized, but 
include the partially ionized species C V-VII, O V-ix, Mg ix-xiii. Si xi-xv, and Fe xvi-xxvii. 

For a given F, the temperature and ionization state of the surface of the slab is expected to 
depend mainly on the value of the ionization parameter, 

(1) 

where F is the total illuminating flux (from 0.01-100 keV) and Hh is the hydrogen number density 
(Ross et al. 1999). 



3. Application to Gamma-Ray Bursts 

If the GRB explodes within a young supernova (SN) remnant, the X-rays from the afterglow 
will illuminate Fe-rich ('^ ten times solar Fe abundance) material, leading to recombination line 
emission by reflection (Vietri et al. 2000). Reasonable values of ^ expected in this scenario are 

^SN ~ 10%8^^r6«^10, (2) 

where L is the source luminosity in units of erg s"^ Hh is the hydrogen density in units of cm~^, 
d is the distance from the burst explosion to the material in units of cm. We adopt the convention 
Q= l(yQx for expressing the physical parameters. The time delay of a few days observed in 
GRB 991216 yields an incident angle for the radiation of i? ~ 45° (Vietri et al. 2000). 
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Altematively, line emission can be produced when a post-burst outflow, possibly magnetically 

dominated, impacts on the walls of a funnel excavated in the stellar envelope (SE) at distances less 
than a light-hour (Rees & Meszaros 2000). Luminosities as high as L ~ 10"^^ erg s"^ are expected 
1 day after the burst, if they are due either to the spinning down millisecond pulsar or to a highly 
magnetized torus around a black-hole (Rees & Meszaros 2000). The ionization parameter in the 
stellar envelope case is 

6E~/510%7^i?rf«ff,i7, (3) 

where /5 < 1 is the ratio of ionizing to MHD luminosity (Rees & Meszaros 2000). The incident 
flux would be deflected along the funnel walls with different incident angles before escaping the 
funnel: for simplicity we assume ~ 45°. 

For these uniform-density slabs we vary ^ by changing the total illuminating flux while keep- 
ing the hydrogen number density hh, the distance from the burst d and the incident angle d fixed. 
Fig. 1 shows the results for illumination by a r=2 spectrum for both of the scenarios described 
above. The illuminating and reflected spectra are displayed as EFe, where Fe is the spectral energy 
flux and E is the photon energy. The models with the highest ionization parameter (^ > 10"^^) are 
excellent reflectors, and show almost no Fe spectral features. This is because the surface layer is 
almost fully ionized, and Fe xxvi does not become dominant until ry > 8. The temperature at 
the surface of the slab is ~ 1.9 xlO^ K, the Compton temperature for the incident spectrum. The 
spectrum drops at ~ 50 keV because we included a sharp energy cut-off in the incident spectrum at 
100 keV^. Compton reflection produces a slight steeping in the reflected spectrum, which mimics 
a power law with F > 2 in the 3-30 keV band. For example, we find r=2.3 in the SN case and 
r=2.14 in the SE case. 

When the illuminating flux is reduced so that 10^^ < ^ < 10"^, Fe Ka emission and Fe K-shell 
absorption features begin to appear. Most of the Ka photons originate at ^ 1 and emerge as 
broad Comptonized lines with weak cores. The H-like Ka photons at 6.97 keV are generated close 
to the surface, but are resonantly trapped and removed from the narrow line core by Compton 
scattering. The He-like inter-combination line at 6.7 keV is not subject to resonant trapping, but 
is generated at such high tt (see panel (b) in Fig. 2) that it is multiply Compton scattered on 
leaving the slab. The broad Comptonized emission features blended into the Compton smeared 
K-shell absorption edge, are an important signature of ionized reflection. The effect of increasing 
the iron abundance at a fixed ionization parameter is illustrated by the dotted line in Fig. 1. The 
Comptonized line and the absorption feature are increasingly significant for an Fe-rich medium. 

With ^ ~ 10^, Fe XXV becomes dominant at tj- 1 (see panel (a) in Fig. 2). Narrow emission 



^For a constant density slab and the values of V considered here, extending the illuminating spectrum to higher 
energies (< 5 1 1 keV) would have Uttle effect on the iron ionization and spectral features (e.g., Ross et al. 1999). 
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line resulting from Fe XXV can now be seen superimposed on the Compton-broadened emission 
bump. The tiny emission feature just above 8.8 keV results from radiative recombination directly 
to the ground level of Fe XX v (see the L46 SE case in Fig 2). In the model with ^se = 10^ the 
Fe emission is suppressed because Fe xvil-xxil dominates close to the surface, but their Ka 
photons are destroyed by the Auger effect during resonance trapping (Ross, Fabian & Brandt 
1996). Finally, for ^se = 10 the reflection spectrum is similar to that of a cold, neutral slab, and so 
the narrow emission line at 6.4 keV is dominant. 

The effect of increasing the incident angle 1?, which is a key ingredient in the pre-ejection 
models because it is responsible for the observed time-delay, is also illustrated in Fig. 1 by the 
dot-dashed line. Radiation that illuminates the atmosphere more directly ionizes more deeply into 
the slab than radiation at grazing incidence. At high flux levels the emergent line features are 
little changed, however, the K-shell absorption features become more prominent for radiation that 
impacts closer to the normal of the slab. 

The ionization structure in the outer layers of the illuminated slab also depends on the incident 
radiation spectrum. Fig. 3 shows the Fe Ka equivalent width (EW) as a function of incident 
luminosity for a series of reflection spectra with F =1.6, 2.0 & 2.4, assuming = 45°. The EWs 
were calculated with respect to the reflection spectrum, and the integration was carried out between 
5.7 and 7.1 keV. As seen in Fig. 1, the EW of the Fe Ka line decreases with L in both GRB 
scenarios. This behavior continues until L is sufficiently low for the narrow Fe Ka emission line to 
be suppressed by the Auger effect. The EWs are generally larger when the illuminating spectrum 
is steeper (i.e., F is greater). The weaker ionizing power of steep spectra allow line emission to 
persist at the highest luminosities, although the EWs still end up quite low when L = L48. In a GRB, 
softer spectra may be important if both the line emission arises from the impact of the continuous 
energy output with the compact remnant and much of the observed continuum emission (with a 
flatter F) comes from the afterglow emission directly. 



4. Discussion 

The emission feature observed ~ 1.5 days after the GRB 991216 burst (Piro et al. 2000) 
had a line luminosity of Lune = 4 x 10"^ erg s"^ and an equivalent width (EW) of ~ 0.5 keV. The 
continuum flux from GRB 991216, measured to have F = 2.2 ± 0.2, in the 1-10 keV band was 
50-100 times stronger than the flux in the line. As is clear from the above discussion, the emission 
feature can be explained by reflection if identified with the recombination Ka line from He-like 
iron at 6.7 keV. 

Taking F Ri 2 in the SE scenario, a continuous ionizing luminosity of ~ 10"*^- 10"^^ erg s"^ (or a 
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continuous wind luminosity of 10''^^ erg s"^ with /3~0. 01-0.1; see Rees & Meszaros 2000) would 
be sufficient to produce the observed line^. However, it is possible that the reflected and incident 
spectra are observed together. With the above parameters, the line luminosity and EW calculated 
for the total spectrum (reflected + incident) are L^ne ~ 2 x IC^"^ erg s"' and EW ~ 0.25-0.7 keV. If 
the illuminating spectrum is a steeper power-law (F > 2), then a smaller fraction of illuminating 
photons lie in the 9-20 keV range which dominates the ionization of Fe, and thus an increase in 
either L, /3 or the Fe abundance is required in order to reproduce the same line strength (Fig. 3). 

Alternatively, an X-ray afterglow with F ~ 2 illuminating the walls of a supernova remnant 
with a luminosity of ~ 10"*^-^ - 10"^^ would produce a line signal with Lune ~ 2 x lO**"* erg s"^ and 
EW ~ 0.2-0.8 keV. However, an essential assumption contained in this model is that the material 
responsible for emission-line features is illuminated by the early afterglow radiation (or the GRB 
itself) with L > L48. Such high luminosities would cause the spectral features resulting from iron 
to disappear or, at best, to be extremely weak (even for larger values of F, see Fig. 3). Moreover, 
this early incident radiation can be harder than a F ~ 2 spectrum with a significant fraction of the 
energy above the 77 — > e^ formation energy threshold, and a high compactness parameter. This 
will cause new pairs to be formed in the originally optically thick scattering medium, an effect 
which amplifies the density of scattering charges and increases the temperature of the illuminating 
material. When pairs are produced in sufficient numbers, the iron Kct emission is suppressed due 
to the decrease in the number of recombinations. These effects will be investigated elsewhere (Laz- 
zati, Ramirez-Ruiz & Rees 2001, in preparation). These problems may be overcome in particular 
source geometries for which lower luminosities and softer spectra are expected at the edges of the 
relativistic outflow. 

The results presented in Fig. 1 show that if the observed Fe Ka line from GRB 991216 is 
identified solely with the H-like line at 6.97 keV, as suggested by Piro et al. (2000), then this is 
inconsistent with emission from a photoionized optically thick slab. Their detailed analysis of the 
Chandra ACIS-S spectrum shows marginal evidence (at the 2.1a level) of an emission feature at 
4.4 ± 0.5 keV. Identifying this feature with the Fe recombination edge with rest energy of 9.28 
keV gives z= 1.1 ±0.1, consistent with the redshift of a H-like Fe Ka line. Nonetheless, it is 
important to emphasize that the observations presented by Piro et al. (2000) do not rule out the 
presence of a He-like Fe feature (or a blend between H and He-like features as shown in Fig. 1). 
Some ambiguity also remains in the rest energy of the emission since Doppler blue-shifts of order 
0.1c are expected in both reflection scenarios (Vietri et al. 2000; Rees & Meszaros 200; Meszaros 
& Rees 2001) and may confuse precise measurements of the line energy. Indeed, for an expansion 
velocity of 0.1c, the emission feature at 4.4 ± 0.5 keV could be attributed to the Fe xxv radiative 



^Under this interpretation, the observed Une width can easily be reproduced by Comptonization for an expansion 
velocity below the Umit of 0.1c inferred by Piro et al. (2000). 
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recombination emission just above 8.8 keV rather than to the Fe recombination edge at 9.28 keV. 

Finally, Co produces both a He and H-like emission line: at 7.242 keV and just above 7.526 
keV, respectively. The exact strength of these features, which we do not include in our calculations, 
could be important in both of the scenarios discussed here and may confuse the identification of 
low signal-to-noise spectral features. 

We have presented calculations of X-ray reflection from Thomson-thick slabs for conditions 
which may arise in the immediate environment of a GRB. Comparisons between putative Fe Ka 
lines detected in X-ray afterglows and those predicted by the computations will be useful in dis- 
tinguishing between the various line emission mechanisms. We expect that more sensitive data on 
X-ray afterglow spectral features will impose strong constraints on the nature of GRB progenitors 
and their environments. 

We thank M. J. Rees, R Meszaros, D. Lazzati, A. Blain and the referee for useful comments 
and suggestions. We are particularly grateful to A. C. Fabian and R. R. Ross for very helpful in- 
sights regarding the calculations. ERR acknowledges support from CONACYT, SEP and the ORS 
foundation. DRB thanks the Commonwealth Scholarship and Fellowship Plan and the Natural 
Sciences and Engineering Research Council of Canada for support. 
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Fig. 1 . — X-ray reflection spectra for illumination with a F=2 power law. In both panels the vertical 
lines are at (from left to right) 6.4, 6.7 and 6.97 keV for the three different Fe Ka line energies, and 
9.28 keV which is the energy of the iron recombination edge. Left panel: illumination of a uniform 
supernova remnant with ten times solar Fe abundance, various values of the afterglow luminosity, 
and = 45°. The ionization parameter is given by ^sn = 10^^48, assuming die=l, «h,io=1 (see Eq. 
2). From top to bottom the equivalent widths and the luminosities (integrated between 5.7 and 
7.1 keV) inferred from the X-ray line features are: 3 eV, 37 eV, 94 eV, 268 eV, 943 eV; 2.5L46, 
2.7L45, 9.6L44, 3.6L44, I.4L44. The dot-dashed line illustrates the effect of increasing the incident 
angle from d = 45° to d = 75° at a constant luminosity. Right panel: illumination of a funnel 
excavated in the stellar envelope with solar Fe abundance and various values of the decaying X-ray 
luminosity from the GRB central engine. In this case ^se - 10^^48, assuming dii=l, nH,ii=l, /3=1 
(see Eq. 3) with -d = 45°. From top to bottom the equivalent widths and the luminosities inferred 
from the X-ray line features are: 3 eV, 62 eV, 1.1 keV, 122 eV, 1.6 keV; 2.6L46, 2.7L45, 2.7L44, 5L42, 
5.2L41. The dotted line shows the reflected spectra for illumination with the same luminosity but 
into a stellar envelope that is ten times more abundant in iron. The line luminosity and equivalent 
width for this model are: 243 eV and 3.7L45. 
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Thomson Depth 

Fig. 2. — Fe ion fractions as a function of Thomson depth produced by the illumination of a 
uniform slab with T=2. The same line type denotes the same species of Fe in both plots. Panel (a): 
illumination of a funnel excavated in the stellar envelope with solar Fe abundance and = 10^. 
(b): illumination of a uniform supernova remnant with ten times solar Fe abundance and ^5^^ =10"*. 
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Fig. 3. — Fe Ka equivalent widths (EW) as a function of both the incident luminosity and F for 
the SE & SN scenarios. -& was taken to be 45° for these models. The EWs were computed from the 
calculated reflection spectra (integrated between 5.7 and 7.1 keV). If one uses the total spectrum 
(incident + reflected) then smaller EWs are obtained. For example, the EWs for r=1.6, 2.0 and 
2.4 when L = L46 (in both cases) are 370 eV, 430 eV & 131 eV (SE) and 43 eV, 145 eV & 805 eV 
(SN). 



